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ABSTRACT 



Several lines of evidence suggest that early-type galaxies might have been very dusty during their 
initial evolutionary stages, characterized by intense star formation activity. The radio selection has 
provided, by far, the richest samples of high-z galaxies, most likely of early type, which may therefore 
yield crucial direct information on the period when the bulk of stars were born. In order to investigate 
the role of dust in these phases, we have analyzed recent observations of a number of high-z radio 
galaxies in the framework of a self-consistent scenario for the evolution of early type galaxies. The 
J> ' available data can be fully accounted for by "opaque" models similar to that already used by Mazzei 

. and De Zotti (1994a) to fit the spectrum of the hyperluminous galaxy IRAS F10214 -1-4724. However, 

there is still considerable latitude for models. Good fits can be obtained with galactic ages ranging 
from 0.05 to 2 Gyr. Crucial constraints may be provided by ground-based sub-mm measurements 
' and by observations with the forthcoming Infrared Space Observatory (ISO). 

in ' 
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C^: 1 INTRODUCTION 

2 ■ The evolution with cosmic time of stellar populations in galaxies can be directly investigated by 

I looking at the spectral energy distribution (SED) of galaxies at different redshifts. "Normal" high-z 

Ci ■ galaxies, however, proved to be remarkably elusive: none of the 100 (out of a total of 104) galaxies 

with measured redshift in the sample of Colless et al. (1993), complete to i? = 22.5, has z > 0.7; 
the highest redshift in the complete i? < 24 sample of Cowie et al. (1991) is z = 0.73; the 73% 
^ I complete redshift survey to i? = 24 by Glazebrook et al. (1995) has discovered only one galaxy at 

. z > 1 {z = 1.108); the measured redshifts of near-IR selected galaxies down to -fC = 20 (nearly 100% 

complete at K < 18 and ~ 70% complete at K = 19-20) are ^ 1, with the exception of one object 
at z = 2.35. 

So far, the only effective method for finding high-z galaxies has been optical identifications of 
radio sources (Chambers, Miley Sz van Bruegel 1988; Lilly 1988; Chambers, Miley & van Bruegel 
1990; Rawlings et al. 1990; Windhorst et al. 1991; McCarthy 1993; Lacy et al. 1994; Dey, Spinrad 
& Dickinson 1995; van Ojik et al. 1994). Of course, the interpretation of the observed SED is 
complicated by the need of disentangling the starlight from the emission of the active nucleus, and of 
properly allowing for effects induced by nuclear activity, e.g. star formation triggered by the radio jet 
as it propagates outwards from the nucleus, which may be responsible for the observed close alignment 
of optical, infrared and radio emissions (see McCarthy 1993, for a recent discussion). 

Nevertheless, high-z radio galaxies can provide unique information on early phases of galaxy 
evolution, when intense star formation occurred. As pointed out by several authors, the failure of 
optical searches to detect primeval early-type galaxies (Djorgovski, Thompson &; Smith 1993; De 
Propris et al. 1993; Parkes, Collins & Joseph 1994) might indicate that they were strongly obscured 
by dust. Franceschini et al. (1994) have shown that under the assumption that ellipticals were 
very dusty during the phases of most intense star formation, a consistent picture can be obtained, 
accounting for the scarcity of high-z optically selected galaxies as well as for the deep counts of galaxies 
at far-IR and radio wavelengths. 

The presence of large amounts of dust in high redshift radio galaxies has indeed been indicated by 
a number of recent observations (Eales & Rawlings 1993; Dunlop et al. 1994; van Ojik et al. 1994; Dey 
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et al. 1995; Carilli 1995). The interpretation of the observed SEDs, however, is not straightforward 
since the commonly used population synthesis models do not allow for dust extinction. 

In this paper we exploit the evolutionary population synthesis models for early type galaxies 
worked out by Mazzei, De Zotti & Xu (1994), taking into account in a self consistent way dust 
absorption and re-emission, to analyze the continuum spectra of high-z radio galaxies for which 
enough data are available; a non-thermal contribution is also included. In § 2 we briefly describe the 
models; in § 3 we present the fits to observed continuum spectra; in § 4 we discuss and summarize 
the main results. 

Throughout this paper we adopt Hq = 50kms~^ Mpc~^. 

2. MODELLING THE SED FROM UV TO FAR-IR 

Following Mazzei et al. (1994), the evolution of the synthetic SED has been linked (albeit roughly) 
to chemical evolution, so that the increased metallicity of successive stellar generations is taken 
into account. The galaxy is assumed to be a close system, i.e. both galactic winds and inflow 
of intergalactic gas are neglected. Matteucci (1992) found that, in the presence of a dynamically 
dominant dark matter component distributed like stars, elliptical galaxies develop winds very late; 
and, in general, we expect winds to be important particularly for low mass galaxies, while we will be 
dealing here with large objects. 

We have adopted Schmidt's (1959) parametrization of the star-formation rate (SFR), ip{t): 



where fg = mgas/nib is the fractional mass of gas in the galaxy {nib is the total mass in baryons), 
assumed to be, initially, unity. Current analyses favour n = 1-2 (Buat 1992), but n = 0.5 has also 
been advocated (Madore 1977). Successful models for early type galaxies require tpo > 100 MQyr~^ 
for n > 1 or 'i/'o > 35 Mq yr~^ for n = 0.5. 

As for the initial mass function (IMF), (t>{m), we have adopted a Salpeter (1955) form: 



with niu = 100 Mq and mi ranging from 0.01 Mq to 0.5 Mq. 

Stellar lifetimes are taken into account (i.e. we have not resorted to the instantaneous recycling 
approximation). The gas is assumed to be well mixed and uniformly distributed. The variations with 
galactic age, T, of the fractional gas mass fg{T) and of the gas metallicity Zg{T) are obtained by 
numerically solving the standard equation for chemical evolution (for details, see Mazzei, Xu & De 
Zotti 1992). 

The synthetic spectrum of stellar populations as a function of the galactic age has been obtained 
using the theoretical isochrones given by Bertelli et al. (1990), as extended by Mazzei (1988) up to 
100 M0 and to an age of 10^ yr. Full details can be found in Mazzei et al. (1992, 1994). 

The effect of dust has been taken into account based on the following assumptions: the dust to 
gas ratio is proportional to the gas metallicity Zg{T)\ the "standard" grain model (Mathis, Rumpl & 
Nordsiek 1977; Draine & Lee 1984), incorporating a power law grain size distribution, holds at any 
time; the size distribution of PAH molecules obeys the power law proposed by Puget et al. (1985), 
while their absorption cross sections and emission properties are those given by Puget &: Leger (1989). 

Corrections for internal extinction have been computed assuming the elliptical galaxies to be 
spherically symmetric, with a density profile described by King's (1966) formula, and adopting the 
interstellar extinction curve given by Seaton (1979) for A < 0.37/xm and by Rieke & Lebofsky (1985) 
at longer wavelengths. The optical depth has been normalized requiring that the far-IR to blue 
luminosity ratio at the present time be equal to the mean value, Lfir/Lb — 7 x 10"'^, derived by 
Mazzei k. De Zotti (1994b) for a complete sample of local elliptical galaxies. Its evolution follows 
directly from that of the gas fraction and of the metallicity. 
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The diffuse dust emission spectrum is modelled taking into account two components: warm dust, 
located in regions of high radiation intensity (e.g., in the neighborhood of OB clusters) and cold dust, 
heated by the general interstellar radiation field. Their relative normalization at the present time 
is obtained fitting the resulting spectrum to the mean local far-IR colours of ellipticals (Mazzei k. 
De Zotti 1994b). The warm to cold dust emission ratio is assumed to be proportional to the star 
formation rate; warm dust is therefore dominant in the early evolutionary phases considered here. 
PAH molecules are mixed with both these components, as discussed in Mazzei et al. (1992). 

The model also allows for emission from hot circumstellar dust, emitting at mid-IR wavelengths. 
Based on the analysis by Ghosh et al. (1986), we assume that the dominant contribution comes from 
stars in the final stage of evolution along the asymptotic giant branch. 

Different choices for the initial SFR, ■0o, and for the index n [eq. (1)], all consistent with the 
observed properties of early type galaxies, result in very different histories for the effective optical 
depth, Teff, defined by LFiR/-Lboi = 1 — exp(— Tes), -Z^boi being the bolometric luminosity of the galaxy. 
This is because Teff is proportional to the product of fg{T) and Zg{T) (as a consequence of the 
assumption of a dust to gas ratio proportional to the gas mctallicity), the first of which is rapidly 
decreasing during the early evolution of the galaxy, while the other is rapidly increasing. 

Therefore, relatively small differences in the evolution of the synthetic spectrum of stellar 
populations, may translate into drastically different amounts of dust absorption and re-emission. 
As illustrated by Fig. 1, models fitting the observed local properties of early type galaxies and 
yielding, at the present time, an effective optical depth equal to the mean value found by Mazzei 
& De Zotti (1994b), may either be optically thin throughout the entire evolutionary history of the 
galaxy (although dust absorption is in any case important during early phases), or become optically 
thick for a short time, or even experience an extended, extremely opaque phase, when essentially 
all light comes out at far-IR wavelengths. Models of the latter type provide a very good fit to the 
observed SED of the hyperluminous galaxy IRAS F10214 -1-4724 (Mazzei k De Zotti 1994a), although 
the recent observational evidences that this source is gravitationally lensed (Matthews et al. 1994; 
Broadhurst & Lehar 1995; Graham & Liu 1995; Trentham 1995; Eisenhardt et al. 1995; Serjeant et 
al. 1995; Downes et al. 1995) complicate the interpretation of the spectral energy distribution. In 
fact, differential magnification of different components having different sizes may change the source 
colours. As pointed out by Downes et al. (1995), it is likely that the near- and mid-IR emissions are 
more amplified than the far-IR emission. 

In the following discussion we will refer only to models with n = 0.5 and n = 1. We stress that, 
in the present data situation, models are necessarily rough and are to be taken only as descriptive of 
a possible astrophysical situation. For example, we regard n = 0.5 models as a simple way to describe 
a situation whereby galaxies become strongly obscured during some evolutionary phase; the physical 
reason for that, however, might be related not to this particular form of the SFR, but, e.g., to merging 
of galaxies with a metal enriched interstellar medium. 

The models successfully match the observed SEDs of local galaxies over four decades in frequency 
(Mazzei et al. 1992, 1994; Mazzei, Curir & Bonoh 1995). As illustrated by Fig. 2a, where 
predictions corresponding to initial SFR's, i/jq, decreasing from values typical of ellipticals (top) to 
values appropriate to irregulars (bottom) are shown, they are also nicely consistent with the observed 
colour distribution of galaxies in the complete K < 20 sample of Songaila et al. (1994); the reddest 
observed colours are close to predictions of dusty models for early type galaxies (solid lines). 

Incidentally we note that, according to our models, evolutionary effects cannot be neglected for 
z > 1; in fact, the colours predicted by models are bluer than those expected by simply redshifting 
the SEDs of template local galaxies, due to the higher SFR. Therefore photometric redshift estimates 
hinging upon the SEDs of unevolved galaxies may be misleading. 

The predicted colour distribution (Fig. 2) is rather sensitive to the adopted cosmological model. 
For given formation redshift Zf, low density models imply considerably redder colours (cf. Fig 2c,f); 
for example, for Zf = 10 [tf = 0.4 Gyr if qo = 0.5 or 1.8 Gyr if qo = 0), we expect values of (B — K) 
and (I - K) of up to ~ 9.6 and 7.7 respectively if q^ = 0, while {B - K) < 8.8 and {I ~ K) < 5.5 
if Qo = 0.5. We note, however, that the colour distribution is also sensitive to Zf, lower values of Zf 
entail bluer (B - K) and (I - K) (Fig. 2b,e). 
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As illustrated by Fig. 3, our models satisfactorily account for the increasing spread of K-band 
magnitudes of radiogalaxies at ^ > 1 in terms of the combined effect of evolution of stellar populations 
and of dust extinction. 

All in all, rather general arguments lead to the conclusion that dust extinction plays an important 
role during the early evolution of elliptical galaxies; there are also indications that a large fraction 
of ellipticals may go through an opaque phase (Franceschini et al. 1994). As pointed out by several 
authors, the inferred ages of high-z radiogalaxies might be reduced if there is significant reddening (cf. 
McCarthy 1993), thus alleviating the problem of galaxy ages uncomfortably high in comparison with 
the age of the universe as well as the difficulty to understand the alignment between optical/infrared 
continua and the (presumably short lived) radio jets. The possibility of reproducing the observed 
SED of high-^ radio galaxies in terms of a young reddened starburst has been questioned, however 
(cf. Chambers et al. 1990). 

In the following section we will analyze the observations of the most extensively studied high 
redshift radiogalaxies to show that, in fact, they are fully consistent with the presence of strong dust 
extinction. 

3. FITTING THE OBSERVED SEDs OF UIGH-z RADIO GALAXIES 

In interpreting the photometric data on radio galaxies it is necessary to take into account that the 
observed luminosity is a mix of stellar and nuclear emissions. As for the shape of the continuum 
spectrum of the latter, the mean SED of local ACNs, derived by Granato & Danese (1994), has been 
adopted. It is also necessary to take into account that low resolution continuum spectra may be 
distorted by contributions of intense emission lines in the broad-band filters (Eisenhardt & Dickinson 
1992; Eales et al. 1993b; Eales & Rawlings 1993). 

We discuss here the four best observed z > 2 radiogalaxies, namely: 4C41.17, the second highest 
redshift galaxy known, B2 0902 + 343, the first z> 3 galaxy discovered, 6C1232 + 39, and 53VK002. 

In general, after subtraction of the emission line contributions to broad-band fiuxcs, the SED 
turns out to be consistent with a young starburst, of age in the range 0.05-0.3 Gyr, assuming that 
the UV continuum is dominated by stellar emission. If, on the other hand, the ACN provides most 
of the rest-frame emission at lOOOA, stellar populations must be older and galactic ages as large as 
1-2 Gyr are derived. 

In any case, there is considerable latitude for models, implying, among other things, that little 
can be inferred from optical/near-IR photometry alone on the amount of dust extinction. Millimeter 
and sub-mm observations of 4C41.17, B2 0902 + 343, and 53VF002 imply large masses of dust and 
the upper limits on 6C1232 + 39 are still compatible with a similar situation. If so, the total baryonic 
masses and the bolometric luminosities must be large: rUb > 10^^ M©, Lboi > 10^'^ L©. 

3.1. 4C41.17 

The published data on this z = 3.8 radio galaxy, the second most distant galaxy known, are shown 
in Fig. 4. The 2.2 /xm fiux is from the line-free Ks magnitude in a 4" diameter beam measured by 
Graham et al. (1994). No correction for contamination by emission lines is possible for I-band and 
R-band fiuxes measured by Chambers et al. (1990); spectroscopic data of the same authors, however, 
indicate that the observed magnitudes can safely be taken as representative of the continuum emission. 
The two points at the shortest wavelengths are from a low dispersion spectrum by Chambers et al. 
(1990). 

This radiogalaxy has been clearly detected at mm and sub-mm wavelengths (Dunlop et al. 1994; 
Chini &; Kriigel 1994). As discussed in both these papers, it is highly implausible that the observed 
emission is of non-thermal origin. The present models straightforwardly account for it in terms of dust 
re-radiation, dominated by warm dust, heated by the hot, young stars. The available data, however, 
do not strongly constrain the warm dust temperature, T^. Shown in Fig. 4 are predictions assuming 
= 60K (as in the case of IRAS F10214 + 4724, cf. Mazzci & De Zotti 1994a) and Td = 46 K. 

Although sub-mm observations were made using a beam size of 16.5" (Dunlop et al. 1994) or of 
11" (Chini &; Kriigel 1994), much larger than the apertures for optical/IR photometry, no aperture 
corrections have been applied to any data. From an image of the region of Lya obscuration, Dunlop 
et al. (1994) conclude that the dust is concentrated in a region of a few arcsec in diameter centered 
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on the flat-spectrum core of the radiogalaxy; any aperture effects should therefore be smaU compared 
with the other uncertainties. 

The 3cr IRAS upper hmits at 12 and 25 //m derived by Dunlop et al. (1994) are at or above the 
upper border of the figure. 

Panel a) of Fig. 4 illustrates fits by "young" models. The galactic ages range from 0.05 to 0.1 
Gyr; thus, 4C41.17 may really be a protogalaxy, in the sense that only a small fraction of stars have 
yet formed; in fact, our models yield a gas fraction fg ~ 0.9. The best fit non-thermal contribution 
(solid line) is ~ 35%-50% of the continuum emission at A = 1000 A. 

The gas metallicity, Zg, and the total baryonic mass, nih, also depend on the lower mass limit, 
mi, of the Salpeter IMF. For 0.01 < mi/MQ < 0.5 we have .0015 < Zg < 0.01, the lowest values 
corresponding to the lowest mi and to the youngest ages, and 1.85 < mb/lO-^^ < 6.5; of course, 
much lower baryon masses be obtained adopting a substantially higher mi. The bolometric luminosity 
is in the range 1.2-3.6 x 10^^ L©. Most of the uncertainty follows from the lack of far-IR data. 
Measurements with the forthcoming Infrared Space Observatory (ISO) should be able both to assess 
the bolometric luminosity and to strongly constrain the dust temperature distribution. 

The internal reddening, in the observer's frame, implied by the models is E{R—K) = 1.0-1.6 (the 
largest value corresponding to the highest value of T), consistent with Eales and Rawlings's (1993) 
estimate, E{R — K) = 1.6, based on line ratios. Models imply a very high SFR, ranging from about 
5000 up to ~ 8OOOM0yr~^, where the lowest value correspond to the highest mi. 

Figure 4b illustrates the predicted SED in the case of a galactic age of 2 Gyr. In this case, the AGN 
must dominate the lOOOA rest-frame emission. The inferred baryonic masses (4 < mb/10^^ M© < 8.6) 
are larger than for "young" models. The residual gas fraction is small {fg ~ 0.0015-0.06) and the 
predicted colour excess, E{R — K), ranges from 0.44 to 0.9 mag with a SFR ranging from about 400 
to 2100 Moyr-^ 

3.2. 55 0902 + 343 

Shown in Fig. 5 are the optical observations of Lilly (1988), the K-band flux, after subtraction of the 
line emission, estimated by Eales et al. (1993b; see also Eisenhardt & Dickinson 1992), the 800 ;um 
upper limit by Hughes, Dunlop & Rawlings (1995) and the detection at 1.3 mm by Chini & Kriigel 
(1994). 

Evidences of strong extinction have been reported by Eales & Rawlings (1993), who estimated 
E{R — K) = 2.1 mag from line ratios, and by Eisenhardt & Dickinson (1992), who pointed out that 
the large [OIIIJ/Lya ratio, the discrepancy between the positions of the radio core and the UV/optical 
emission peaks, and the strong spatial variations in the R — K colour could all be interpreted as due 
to internal extinction. 

Panels a) and b) of Fig. 5 show fits by "young" and somewhat older models, respectively. Fits 
were found only for ages of less than 1 Gyr. 

3.3. ^C1232 + 39 

Deep optical, infrared, and radio images of this z = 3.22 radio-galaxy have been obtained by Eales 
et al. (1993a). The K-band continuum flux is very uncertain: by comparing their spectroscopy with 
published photometry, Eales & Rawlings (1993) derived an estimate of 34 ± 10 //Jy, whereas there is 
no evidence of a continuum flux in their infrared spectrum. Plotted in Fig. 6 is the 3(T upper limit 
(37 //Jy) on the continuum emission measured by Eales & Rawlings (1993) from the spectrum. Eales 
et al. (1993a) point out that the intensity peaks of B and R images are offset from the K-band peak, 
which coincides with the center of the radio image; they suggest that this difference in structure may 
be due to dust obscuring the emission from the center of the galaxy at the shorter wavelengths. The 
upper limits at 800 //m and at 1.3mm are from Hughes et al. (1995) and Chini &; Kriigel (1994), 
respectively. 

No estimates of the non-thermal contribution to the observed fluxes are available. Again, "young" 
models, suggesting 0.05 < r(Gyr) < 0.2, imply a minor (~ 10%) contribution from the AGN at 
Arestframe = 1000 A; if this Contribution is increased to, say, ~ 60%, fits require a somewhat older age, 
but stiU < 1 Gyr, as in the case of B2 0902 + 343. 

Because of the lack of a reliable estimate of both the dust extinction and the far-IR emission, 
models allow the possibility of either small or large (up to E[R — K) = 1.7 mag) dust extinction. 
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This z = 2.39 radio-galaxy (Fig. 7) is the most distant so far discovered by the LBDS survey 
(Windhorst et al. 1991). HST WFC V-band and I-band images (Windhorst et al. 1992; 1994) show 
a radial intensity profile consistent with that of an elliptical galaxy, plus an unresolved central source 
(AGN) containing 30 ±10% of the total light. Under the plausible assumption that the galaxy surface 
brightness is monotonic with radius, Windhorst et al. (1992) conclude that the AGN contributes 
at most 20% in V (A^est — 1600 A) and 25% in I (Arest — 2600 A), with a secure upper limit of 
37%. The ground based spectroscopy of Windhorst et al. (1991) indicates a non thermal fraction at 
Arest = 1450 A of 35 ±15%. 

The near-IR spectrum obtained by Eales & Rawlings (1993) is very noisy and the correction 
for the line contribution in the K-band is correspondingly very uncertain. We have adopted the 
line-corrected magnitude given in Table 2 of Windhorst et al. (1994). 

Assuming that the observed low value of the ratio Lya/Ha is due to dust extinction, Eales &; 
Rawlings (1993) obtained a quite substantial dust reddening, E{R — K) ~ 2.4. The presence of a large 
amount of dust may be supported by the (marginal) detection of this galaxy at 800 fiva by Hughes et 
al. (1994). Our models are consistent with 0.75 < E{R ~ K) < 1.7. Good fits to the observed SEDs 
are obtained for ages down to 0.1 Gyr. However, its a^/^-like light profile observed by Windhorst et 
al. (1994) indicates dynamical relaxation and, therefore, this galaxy cannot be too young. In Fig. 
7 we show models with ages of 0.6 and 1 Gyr. The assumed AGN contribution (~ 40% in V) is 
consistent with the results of ground based spectroscopy by Windhorst et al. (1991) and close to the 
upper limit derived from HST images by Windhorst et al. (1992). Note, however, that the quality of 
the fit is not crucially dependent on the AGN emission. 

4. CONCLUSIONS 

Our analysis of the observed SEDs of the best studied z > 2 radiogalaxies shows that all the 
ingredients (stellar population synthesis, dust extinction, contribution of the AGN) need to be taken 
simultaneously into account. Ignoring one of them may lead one astray. 

In general, the available data are consistent with galactic ages ranging from ~ 0.1 to 2 Gyr. 
Young models require small AGN contributions in the UV; if, on the other hand, stellar populations 
are 1-2 Gyr old, the UV light may be mostly produced by the AGN. 

Young models are consistent with (although do not require) large amounts of dust extinction 
such as that indicated by the 800 //m flux detected from 4C41.17 (Dunlop et al. 1994), by mm 
measurements by Chini & Kriigel (1994), and by the possible detection at 800 /im of 53W002 by 
(Hughes et al. 1995). The presence of substantial amounts of dust is also suggested by the low values 
of Lya/Ha and Lya/H/J ratios measured for several high-z radio-galaxies (Eales & Rawlings 1993). 

The existence of a dusty phase in the evolution of early type galaxies has been advocated to 
explain the remarkable scarcity of high-z galaxies in faint optically selected samples as well as the 
strong cosmological evolution in the far-IR suggested by deep IRAS counts at 60 jim. (Pranceschini et 
al. 1994). 

High sensitivity sub- mm and mm observations such as those carried out by Dunlop et al. (1994), 
Hughes et al. (1995) and Chini k. Kriigel (1994), as well as measurements with the long wavelength 
photometer (PHOT) on the forthcoming ISO satellite might detect the corresponding dust re-emission 
or, at least, set significant constraints on the models. 
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Figure captions 



Fig. 1. Evolution of the effective optical depth of dust (defined in § 2), of the gas metallicity and of 
the gas fraction for various choices of the parameters controlling the SFR and of the lower mass limit 
of Salpeter (1955) IMF. 

Fig. 2. Distribution of {B — K) and (/ — K) colours as a function of redshift for galaxies in 
the complete {K < 20) sample of Songaila et al. (1994); stars are for galaxies lacking redshift 
measurements, arbitrarily set at z = 5. The lines show the relations implied by the present models 
for different choices of the SFR, of the redshift of galaxy formation Zf, and of qo- In panels a and d, 
qo = 0.5 and zj = 10 (except for the model represented by the dot-dashed line, which has Zf = 2.5); 
in panels b and e, qo = 0.5 and Zf = 5 (but the dot-dashed line has Zf = 1); panels c and / are for 
qo = and = 10 (but the dot-dashed line has Zf = 2.5). In each panel, the solid line refers to 
ipo = lOOM0yr~^ and n = 0.5, i.e. to dusty early-type models; the dashed line, to ipo = 5OM0yr~^ 
and n = 0.5; the dotted hnc to V'o = 100 Mq yr~^ and n = 1; the dot-dashed line to ■00 = 3.5 Mq jr~^ 
and n = 1, appropriate for late type galaxies. 

Fig. 3. Predicted vs observed K-z diagram of radiogalaxies for different values of the galaxy 
formation redshift, Zf, and of qo. In each panel, the solid line refers to ^po = lOOMQyr"^ and n = 0.5 
(dusty early-type galaxies); the dashed line, to ipo = 5OM0yr~^ and n = 0.5; the dotted line to 
ijjQ = lOOMQyr"-'^ and n = 1. The data points are those collected by Bales et al. (1993a; their Fig. 
3). The open symbols correspond to line-corrected magnitudes. 

Fig. 4. Spectral energy distribution (SED) of 4C41.17. All models have n = 0.5 [cf. eq. (1)]. The 

solid line shows the adopted non-thermal SED (see text). In panel a, model A has age T = 0.1 Gyr, 
dust temperature = 46 K, total baryonic mass mi, = 6.5 x 10^^ Mq, bolometric luminosity Lboi = 
1.2 X 10^^ Lq, E{R — K) = 1.03, gas fraction fg = 0.90, gas metallicity Zg = 1.5 x 10~^, non-thermal 
fraction at lOOOA (rest-frame) /nt = 0.45; model B has T = 0.1 Gyr, Ta = 46 K, nib = 8.4 x 10^^ 
Lboi = 3.2 X 10^3 Lq, E{R - K) = 1.3, fg = 0.91, Zg = 3.6 x lO'^, /nt = 0.42; model C has 
T = 0.05 Gyr, Td = 60 K, rUb = 1.9 x 10^^^^, Lboi = 3.6 x IO^^Lq, E{R - K) = 1.63, fg = 0.88, 
Zg = 0.01, /nt = 0.35. All models in panel h have T = 2 Gyr and /nt = 1- The values of ■00 and 
mi are the same as in panel a. Model A has = 60 K, mb = 8.6 x 10^^ M©, Lboi = 2.8 x 10^^ L©, 
E{R -K) = 0.90, fg = 0.06; model B has Ta = 60 K, rub = 6.4 x lO^^^^, Lboi = 2.1 x lO^^L©, 
E{R - K) = 0.90, fg = 0.06; model C has Ta = 46K, = 4 x IO^^Mq, Lboi = 2.45 x lO^^L©, 
e{r - K) = 0.44, fg = 0.0015. The data points are described in the text. 

Fig. 5. SED of B2 0902 -|- 343. Ah models have n = 0.5. In panel a is = 46K for all models. 
Model A has T = 0.3 Gyr, = 5.3 x IQ^'^Mq, Lboi = 2.3 x IO^^Lq, ^{R - K) = 1.30, fg = 0.77, 
Zg = 1.3 X 10-2, /nt = 0.05; model B has T = 0.1 Gyr, nib = 2.5 x lO^^ ^^^^ ^ o.97 x lO^^ 
E{R-K) = 1.32, fg = 0.91, Zg = 4x10-3, /nt = 0.1; model C has T = 0.05 Gyr, mb = 1.5x10^2^^^ 
Lboi = 2.9 X 10^3 Lq, EiR -K) = 1.07, fg = 0.88, Zg = 0.01, /nt = 0.1. 

In panel 6, /nt = 0.6, ipo = lOOMQyr ^, and mi = 0.1, = 60K for all models. Model A has 
T = 0.6 Gyr, mb = 7.9 x lO^^^^^ ^.boi = 3.5 x lO^^L©, E{R - K) = 1.4, fg = 0.76, Zg = 0.028; 
model B has T = 0.6 Gyr, = 6.4 x lO^^M©, Lboi = 2.9 x lO^^L©, E{R - K) = 1.3, fg = 0.76, 
Zg = 0.028. The data points are described in the text. 

Fig. 6. SED of 6C1232 + 39. Models A and B have n = 0.5, model C has n = 1. In panel a 
is Td = 60 K for all models. Model A has T = 0.1 Gyr, mb = 3 x IO^^^q, Lboi = 1-3 x IO^^Lq, 
E{R-K) = 1.5, fg = 0.9L Zg = 4x10-^, /nt = 0.13; model B has T = 0.05 Gyr, mb = 2x10^2^0, 
Lboi = 4.2 x IO^-'^Lq, E{R ~'K) = 1.7, fg = 0.88, Zg = 0.01, /nt = 0.1; model C has T = 0.2 Gyr, 
mfo = 2 X 10" M0, Lboi = 2.6 x IO^^l®, E{R - K) = 0.24, fg = 0.20, Zg = 0.08, /nt = 0.09. In 
panel 6, models A and B have /nt = 0.6 and Td = 60 K. Model A has T = 0.8 Gyr, mb = 4:X 10^^ Mq, 
Lboi = 2 X IQi'^Lo, E{R - K) = 1.35, fg = 0.09, Zg = 0.18; model B has T = 0.4 Gyr, mb = 
4 X 10" M©, Lboi = 3 X 10^2 L©, E{R - K) = 0.2, fg = 0.05, Zg = 0.15; model C has T = 0.6 Gyr, 
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nib = 2.2 X 10^2 Lboi = 1-0 x 10^^ ^{R - K) = 0.75, fg = 0.57, Zg = 0.03. The data points 
are described in the text. 

Fig. 7. SED of 53VF002. AU models have n = 0.5 and = 100. The sohd Unc shows the adopted 
non-thermal contribution. Model A has T = 1 Gyr, rrih = 1.9 x 10^^ M©, Lboi = 3.4 x 10^^ L©, 
E{R - K) = 1.35, fg = 0.30, Zg = 0.024; model B has T = 0.6 Gyr, = 1.2 x IO^^m^^ 2>boi = 
5.4 X 10^2 -K) = 1.5, fg = 0.57, Zg = 0.028; model C has T = 1 Gyr, = 1.2 x 10^^ y^^^ 

Lboi = 5.3 X 10^2 ^(^ji -K) = 0.84, fg = 0.36, Zg = 0.051. The data points are described in the 
text. 
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